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Abstract 


Shielding strategies for extended stays in space must include a comprehen- 
sive resolution of the secondary radiation environment inside the spacecraft 
induced by the primary, external radiation. The distribution of absorbed dose 
and dose equivalent is a function of the type, energy and population of these 
secondary products. A systematic verification mid validation effort is under- 
way for HZETRN, which is a space radiation transport code currently used by 
NASA. It performs neutron, proton and heavy ion transport explicitly, but it 
does not take into account the production and transport of mesons, photons 
and leptons. The question naturally arises as to what is the contribution of 
these particles to space radiation. The pion has a production kinetic energy 
threshold of about 280 MeV. The Galactic cosmic ray (GCR) spectra, coinci- 
dentally, reaches flux maxima in the hundreds of MeV range, corresponding to 
the pion production threshold. We present results from the Monte Carlo code 
MCNPX, showing the effect of lepton and meson physics when produced and 
transported explicitly in a GCR environment. 


1 Introduction 

Protecting astronauts from the harmful effects of space radiation is an important pri- 
ority for human space flight [1, 2], The shielding strategies for an extended stay in space 
must include a comprehensive resolution of the internal radiation environment induced by 
the primary, external radiation. Primary radiation particles undergo atomic and nuclear 
interactions as they pass through matter producing secondary particles. The radiation 
downstream consists of modified particle fluxes. The interactions of various radiations 
with matter are unique and determine their depth of penetration in matter. This conse- 
quently impacts the type and amount of shielding needed for radiation protection. Hence, 
for optimal radiation shield design, a complete characterization of the secondary radiation 
products is necessary. 

The space radiation environment is comprised of energetic particles produced from 
three sources, each with a characteristic spectrum. Firstly, solar particle events (SPE) 
consist primarily of protons emitted from the Sun during coronal mass ejections and solar 
flares. These events are rare, but when they occur, they can inflict a potentially lethal 
dose of radiation to astronauts if no protective measures are undertaken. They are also 
of great concern for the stability of electronic devices. Energies often reach hundreds of 
MeV and can even extend into the GeV region. Secondly, Galactic cosmic rays (GCR) 
are protons and heavier nuclei thought to be emitted from supernovae explosions within 
Milky Way galaxy and accelerated to the vicinity of the Solar system. The GCR particles 
have energies from a few MeV per nucleon up to the ZeV region (Zetta eV = 10 21 eV). The 


1 



GCR spectra reaches flux maxima in the hundreds of MeV range and so the MeV - GeV 
region is most important for space radiation. Radiation dose from nuclei is approximately 
proportional to Z 2 , where Z is the ion charge, and so GCR ions from Hydrogen to Nickel 
are of most concern [1], Beyond Nickel, the particle flux is much smaller and is generally 
ignored. Thirdly, the geomagnetically trapped particles are protons and electrons confined 
by the magnetic field of Earth. There are two distinct regions called the inner and outer 
Van Allen radiation belts. Protons and electrons are found in both belts. The proton 
energies range up to 100 MeV and beyond. The electrons energies range from 100 keV 
and beyond. For the purpose of this study, we confine our external radiation environment 
to GCR. Parallel studies are underway for the SPE and trapped belt environments. 

A comprehensive radiation shielding design study requires characterizing the primary 
radiation and the resulting secondary radiation. HZETRN (High Z and Energy TRaNs- 
port, where Z is the charge), is a space radiation transport code currently used by NASA 
to characterize the space radiation environment inside spacecraft for human exposures 
[2, 3, 4], It performs neutron, proton and heavy ion transport explicitly, but it does not 
take into account the production and transport of mesons, photons and leptons [4, 5]. 
The question naturally arises as to how important these particles are with respect to 
space radiation problems. Compared to heavier mesons, pions are the lightest meson and 
are therefore produced more copiously in cosmic ray interactions. Subsequently, through 
decay or other interactions, pions produce hadrons, photons and leptons. 

The primary purpose of this paper is to quantify the difference in dose when pions 
are produced and transported explicitly in a GCR environment. In section 2, we briefly 
summarize the relevant theory of hadron physics. In section 3, the radiation transport 
code MCNPX (Monte Carlo N-Particle extended) and the organization of the input 
parameters for the comparative test are described. The results are reported in section 4. 
The main differences in absorbed dose, and flux due to meson interaction processes are 
discussed separately. Section 5 summarizes the findings and provides some concluding 
remarks. 


2 Hadron physics 

We begin with a review of the relevant physics pertaining to pion production. The 
Standard Model of particle physics [6] describes our universe in terms of fundamental par- 
ticles, called quarks and leptons, which interact via the electromagnetic, strong, or weak 
force. The theories of these forces are quantum electrodynamics (QED), quantum chro- 
modynamics (QCD), and quantum flavordynamics (QFD), mediated by particles called 
photons (7), gluons, and W ± ,Z° bosons respectively. These force mediating particles 
are called guage bosons. Gravity is not included in the Standard Model. There are 6 
quarks called up (w), down (d), strange (s), charm (c), bottom ( b ) and top ( t ). There 
are 6 leptons called electron (e), muon (yu), tau (r), electron neutrino (z/ e ), muon neutrino 
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(u^) and tan neutrino (z/ r ). Antiparticles are denoted with an overline symbol, such as 
the anti-up quark (u). The name for these different types of quarks and leptons is called 
flavor, and weak interactions can produce flavor changes such as the decay of a muon into 
an electron. This is why the theory of weak interactions is called quantum flavor dynamics. 

Quarks and gluons carry another quantum number called color, and hence the name 
of chromodynamics describing the strong interactions, which involve color changes. The 
values of this color charge are called red (R), green (G), blue ( B ), anti- red (R), anti- 
green (G), and anti-blue (B). A red quark can change to a blue quark by emitting a 
gluon carrying red and anti-blue quantum numbers. Color is like magnetism, in that just 
as an isolated North or South magnetic pole cannot exist by itself, so too color cannot 
exist isolated by itself. Therefore, isolated quarks and gluons are not found in nature. 
Particles which do not carry color, such as leptons, photons, and W ± , Z° are found as 
isolated particles. Quarks and gluons can only be found in bound states, such that the 
total color charge is zero, or “white”. The combination of red, green and blue gives 
“white” . Mathematically, standard combinations are R + G + B = R + G + B = 0 and 
also R + R = G + G = B + B = 0. The particles consisting of quark bound states are 
called hadrons, which come in two varieties called baryons and mesons. The baryons, 
such as a proton or neutron, consist of 3 quarks, each with a different color (R,G,B). 
The mesons, such as a pion or kaon, consist of 2 quarks carrying opposite color such as 
(R, R). Actually, hadron wave functions are a linear combination of all allowed values of 
color. For example, a proton is made of the flavors uud and a neutron is made of udd, 
with each of the 3 quarks carrying different color. 

The positive pion (7r + ) is made of ud and the positive kaon (A' + ) is made of us, with 
each quark carrying (a linear combination of) opposite color. The spectrum of hadrons 
consists of all possible allowed combinations of bound states of quarks. These bound 
states occur in both ground and excited states, resulting in a large number of possible 
hadrons with a large variety of masses, determined from the quark masses and the energy 
of the excited hadron states. The lowest mass baryon is the nucleon, which has two charge 
states: namely the proton (p) and neutron (n). The lowest mass meson is the pion, which 
occurs in three different charge states (7r ± , 7r°) . 

Cosmic ray interactions with matter include high energy proton - nucleus and nucleus 
- nucleus collisions; whereby, a nucleus may break up into its constituent pieces, producing 
lighter nuclei in the final state through nuclear fragmentation. There are a large variety of 
different possible reactions; one example is 56 Fe+ 26 Al -^ 55 Fe+ 26 Al + n. The constituents 
of nuclei are protons and neutrons and therefore a nucleus - nucleus collision is basically a 
collection of nucleon - nucleon collisions. Most space radiation studies include the baryons 
but neglect the radiation dose produced by mesons, leptons and their decay products [6] . 
Because the pion is the lowest mass hadron, it is the most copiously produced particle in 
the nucleon - nucleon collisions that occur in cosmic ray nuclear interactions. The heavier 
mesons are produced in fewer numbers. The question naturally arises as to what is the 
contribution of these particles to space radiation. 
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Threshold energies for several pion producing reactions in proton - proton (pp) col- 
lisions are listed in Table 1. The threshold for n° production is at a kinetic energy of 
280 MeV. Double pion production begins at 592 MeV. The GCR spectra, coincidentally, 
reach flux maxima in the hundreds of MeV range, corresponding to the pion production 
threshold. Hence, in a GCR environment (proton peak at about 360 MeV), we expect 
to see both neutral and charged pion production. Further, we also expect to see other 
nucleon - nucleon reactions producing pions, as listed in Table 2. 


Table 1: Kinetic energy thresholds (MeV) for proton - proton (pp) reactions. Particle 
symbols are proton p, neutron n, deuteron d and pion n. 


Final state 

Threshold 

(MeV) 

ppn° 

280 

dn + 

288 

pmr + 

292 

pmr + TC° 

592 

ppn + ii- 

600 

pp7T + 7r~7r 0 

920 


Once produced, the neutral pion will decay immediately; whereas, the charged pions 
will travel some distance before they decay. The pion decay modes and lifetimes (r) are 
listed in Table 3. In Table 3, the lifetime r has been listed as well as the quantity cr. 
Both of these are given in the rest frame of the decaying particle. The ir° decays via the 
electromagnetic interaction, whereas the 7r ± decay via the weak interaction (signaled by 
the presence of a neutrino in the final state), into a muon and neutrino. The produced 
muon subsequently undergoes a weak interaction decay into an electron and neutrinos. 
The electromagnetic interaction, being stronger than the weak interaction at this energy 
scale, accounts for the much shorter lifetime of the 7T° compared to the charged particles. 
In terms of quarks and gnage bosons, the negative pion and muon decays are as follows: 

7 t~ = du — > W~ — > v~ ( 1 ) 

and 

p~ — > e"Z° — > e~V e v ll (2) 
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Table 2: Some nucleon- nucleon reactions producing pions. 


p + p — > p + p 

n + p + 7i + 

7T° + p + p 

^ t + + p + n 

7T + + d 

7T~ + p + p + 7T + 

n + n — > p + n + 

n + n 
7T° + n + n 
7T + + 7i~ + n + n 
7 t~ + p + n 
tt~ + d 

p + n — > p + n 

n + p 

7T° + p + n 
7T° T d 
7T + + n + 71 
7T _ + p + p 


with the positive particles decaying via charge conjugate reactions [6]. The copious pro- 
duction of electromagnetic particles, such as electrons, positrons, photons and muons is 
the main source of electromagnetic (EM) cascades. For example, positrons produced from 
a pion interaction will further interact with electrons to produce photons. Photons can 
subsequently produce e + e _ pairs, resulting in an EM cascade. 

The quantity cr gives an idea of how far the particle will travel before decaying. Of 
course, in the lab frame, which is the target frame or spacecraft wall frame, the lifetime 
will be longer, due to time dilation, and the distance will therefore be longer. Thus, cr 
is actually a minimum distance, but it gives a rough idea of the range of a particle. In 
the flux versus depth and dose versus depth curves, we expect to see the effect of the 
neutral pion decay through the build up of photons. Since the lifetime of charged pions 
is large, before they decay they will interact with other nucleons, producing pions and 
nucleons in the final state. Examples of such reactions are listed in Table 4. These pion 
interactions will impact the overall flux of light ions and baryons. Hence, in a GCR 
transport calculation using a transport code that includes hadron (baryons and mesons) 
and lepton calculations explicitly, we expect to see production of EM particles and EM 
cascades. 
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Table 3: Decay modes [6]. The mean lifetime is given by the symbol r and cr is the speed 
of light multiplied by the mean lifetime. Particle symbols are pion n, muon p, photon 7 
and neutrino u. 


Particle 

Rest Mass 
(MeV/c 2 ) 

Decay 

Mode 

r 

(sec) 

CT 

(m) 

7T° 

135 

77 

8.4 x 10~ 17 

25 x 10“ 9 

7 T + 

140 


2.6 x 10~ 8 

7.8 

7T — 

140 


2.6 x 10' 8 

7.8 


106 

e + u e U^ 

2.2 x 10~ 6 

660 

pT 

106 

e"F e z/ M 

2.2 x 10~ 6 

660 


3 Radiation transport 

The interactions of the incident primary radiation with the target nuclei (shielding or 
human tissue) result in a modified radiation field downstream. A detailed and accurate 
characterization of the radiation field in terms of absorbed dose (D), dose equivalent (H) 
and other end-point quantities, is necessary to assess the short and long term response of 
biological systems and electronics to the exposed radiation. 

The primary radiation particles from a GCR spectrum travel through the spacecraft 
material producing secondary particles. The nuclear fragmentation processes result in a 
cascade of secondary products: leptons, baryons, mesons, photons, electrons, positrons, 
light (A <4) and heavy ions (A > 4). For a comprehensive transport calculation, all 
the secondary products should be modeled explicitly. However, a fully coupled transport 
calculation that includes all possible secondary products is typically not performed due 
to two main reasons. First, these types of calculations traditionally impose considerable 
penalty on time-to-solution for even simple geometries. Second, due to their lower weight- 
ing factors (quality factors), the dose equivalent contribution from leptons and mesons 
is perceived to be negligible when compared with the dose contributions from light ions, 
heavy ions and baryons. Hence, in most transport calculations, certain production and 
transport channels are not explicitly defined, or are turned off to speed up the calculations. 
In many cases, these assumptions have been valid. However, review of the particle theory 
suggests that for the GCR environment, the meson, lepton and photon contributions may 
merit further investigation. 

A systematic verification and validation (V&V) effort, including comparisons to Monte 
Carlo (MG) codes, has begun for HZETRN [7]. Several general purpose MC codes are be- 
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Table 4: Some pion - nucleon reactions. 


7 T° + p 


7T° + n 


7T + + p 


7T + + n 


7T +p 


n + n 


p + 7T° 
n + 7r + 
7T° T p 
7r + + n 


7 T~ +p 

+ 7T+ 

p + n~ 


n + 7T° 


7T° + n 


7 r + + n 

+ 7T^ 

7T _ T 7T 

+ + n 

P + 7T + 


n + 7 r + 

+ 7T + 

7T° T p ■ 

T 7T + 


+ 7T+ + 7 T + 


+ 7T+ 


+ 7T~ 

+ 7T _ 

7T U + n + 7T _ 

7T + + U + 7T _ + 7T 

7l~ + n 


7T + + p 
7 T~ +p 
P + 7T° 
n + 7r + 
7T° T p 

7r + + n 

7 T~ +p 
P + 7T“ 
n + 7T° 
7T° + n 

7r + + n 

7T — T p 
P + TT~ 

n + 7T _ 


ing used for V&V of HZETRN. These include: Particle and Heavy Ion Transport code Sys- 
tem (PHITS), FLUKA, MCNPX, High Energy Transport Code-Human Exploration and 
Development of Space (HETC-HEDS), MARS, GEANT4 and others [8, 9, 10, 11, 12, 13]. 
These MC codes are routinely used in V&V activities of simulated and measured data in 
space environments, and are particularly important to simulate coupled EM and hadronic 
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cascades. Mesons, leptons and photons are currently not included in the HZETRN trans- 
port code [3, 6], however, most general purpose MC codes include meson, lepton, photon 
and other secondary particle physics and transport. Each code adopts approximations 
and simplified treatments to speed up the computation, which can introduce systematic 
effects. There have been numerous investigations that studied the differences in these 
codes and the implications of embedded assumptions to the final results. A number 
of investigators who have studied secondary particle production using MC codes, have 
noted agreements and disagreements in particle fluxes, dose and other related quantities 
[14, 15, 16, 17, 18, 19]. The differences are attributed to the distinction in the assumptions 
made in the transport calculations, or the differences in the high energy physics models 
that different codes use. In our literature review, we did not find an effort that quantifies 
the effects of the meson, lepton and photon transport to space radiation. Some peripheral 
references were found. For example, Aghara et al. [14, 15] quantified photon production 
rates from various proton beam energies impinging into Martian and Lunar regolith tar- 
gets. They noted an increase in the photon production rate at proton kinetic energies E 
> 250 MeV. They showed that two photons are produced for one proton source particle 
(E > 250 MeV). This study did not explicitly link the photon production to pions; how- 
ever, it demonstrated that the rate of photon production increases as the incident energy 
approaches the pion production threshold energy. Similarly, a number of researchers have 
compared simulation results with measurements [20, 21, 22], Most of these studies were 
geared towards code V&V and benchmarking. In summary, the literature review reveals 
that there has not been a conclusive effort to quantify the effects of leptons and mesons 
and their contribution to absorbed dose in a GCR environment. Hence, calculations per- 
formed in this paper are designed specifically to quantify the lepton and meson fluxes and 
their contribution to the total absorbed dose in a GCR environment. 

3.1 MCNPX transport 

The MCNPX code package version 2.6F is used for all the calculations presented in 
this paper. It was chosen for three main reasons; the MCNPX code has been used for 
space radiation research and is well benchmarked, it has an extensive statistical analysis 
package for tally results, and it provides the user with the ability to selectively turn on 
and off specific particle physics. We now discuss some pertinent features of MCNPX. 

MCNPX was developed at Los Alamos National Laboratory, and is a three dimen- 
sional, general purpose Monte Carlo transport code. It was developed by combining 
LAHET nuclear physics models and high energy transport to the MCNP coupled neu- 
tron, photon, electron code [10, 23]. MCNP and MCNPX publicly released versions are 
available through the Radiation Safety Information Computational Center at Oak Ridge 
National Laboratory (RSICC) and the Organization for Economic Cooperation and Devel- 
opment / Nuclear Energy Agency (OECD/NEA). The last MCNPX general distribution 
version, version 2.5.0, only transported particles with atomic weight of four and below. 



The MCNPX code package has been used by a number of investigators to model the 
complex geometry of spacecraft, and extraterrestrial bases (planetary habitats) in recent 
years [14, 15, 17, 22], Although these calculations were limited to baryons and light ions 
up to alpha particles, they provide valuable benchmark verification. The beta version 
2.6F, released in March 2008, added the capability to perform calculations using approx- 
imately 2205 heavy ion particles [24], This new capability in MCNPX enabled a fully 
coupled, GRC spectral ion simulation. These capabilities are available in the new general 
distribution version 2.6.0 released in March 2008 [25]. 

The MCNPX code package uses the MCNP4C transport package for electron, photon 
and neutron transport coupled with physics packages for the lepton, baryon, meson and 
light ion transport, ft includes continuous energy transport of 34 different particle types 
(neutron, photon, electron, 5 leptons, 11 baryons, 11 mesons, and 4 light ions) up to 1 
TeV per nucleon. 

The new heavy ion transport and physics is implemented such that users can continue 
to use the powerful MCNPX capability of modeling a variety of sources (volume, surface, 
distributed, etc.) and tally options [10, 25]. It also includes an extensive tally analysis 
package to determine the statistical error in the simulation results. The statistical analysis 
of the results, a powerful feature in MCNPX, is one of the most critically important 
requirements for high confidence in MC calculations. 

In order to study the effect of pions to the absorbed dose, a code is needed that can 
turn on and off various physics options without changing the problem setup. MCNPX 
provides the user with this capability. For example, when the pion channel is turned off, 
the code still conserves the total energy and charge distribution. It just does not produce 
pions and transport them. The subsequent pion interactions would not be considered 
through the phase space. Hence, if the pion physics and transport is not important to a 
problem, then a statistically significant difference would not surface in the results. In some 
cases, the secondary particle would not show a significant difference. For example, the 
secondary particle may not be produced in large numbers or may have a very low interac- 
tion probability. In both of these cases, it would lead to a negligible difference in the final 
results with the channel on and off. It is important that physics and transport changes 
for a particular particle are imposed without making any changes to the experimental 
setup (i.e. source, geometry phase space). Any difference observed should be entirely due 
to the transport and interaction of the particle in the phase space. Lastly, it is impor- 
tant that the physics models are validated against measured data or compared to another 
code where data is not available. Collot et al. [19] compared pion yield and characteristic 
spectra between FLUKA and MARS with some measurements. They claim reasonable 
agreement between data and code simulation results. MCNPX uses the same Los Alamos 
Quark-Gluon String Model (LAQGSM) and Cascade-Exciton Model (CEM) models used 
in MARS. MCNPX is also well benchmarked for light ions (A<4), neutron, photons, and 
electrons. Hence, the addition of a heavy ion source capability makes MCNPX ideal to 
investigate the impact of the lepton and meson physics for this study. 
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For the heavy ion transport calculations, MCNPX 2.6F uses the new modified version 
of the LAQGSM event generator to compute individual interactions with heavy energetic 
particles and heavy targets [24, 26, 27, 28]. A brief overview of the heavy ion physics 
implementation in MCNPX is provided in the following section. 

3.2 MCNPX heavy ion physics 

MCNPX utilizes nuclear cross section libraries to transport neutrons, photons, pro- 
tons and electrons where measured data is available. The measured cross-section data is 
available for some isotopes up to a kinetic energy of approximately 150 MeV, with the 
availability varying based on isotope. The Los Alamos data group updates the cross- 
section data libraries, and once available, these data tables are included into the MCNPX 
cross section table libraries. For particles and energies where data tables are not available, 
the code uses physics models [10, 29]. The LAHET code system calculates the genera- 
tion (spatial and mass distribution) and point of origin of recoiling particles ejected from 
nuclear reactions, including intranuclear cascades and subsequent fission or evaporation 
[30, 31]. 

The energy threshold for transition from library and model physics is nuclide de- 
pendent, and the transition is handled automatically by a mix-match algorithm in the 
code. Table 5 shows a condensed description of the transition of event generator codes 
using mix-match capability. Details of the mix-match capability can be found in the Los 
Alamos report by Hendricks [29]. The user can override these transitions as needed and 
when possible using physics card options. 

Within the hadron model physics region, MCNPX defines two energy regimes based 
on projectile energy; E < 3.5 GeV and E > 3.5 GeV. The particle interactions in the 
energy region between the library physics and E < 3.5 GeV is handled by five different 
intranuclear cascade (INC) models: Bertini, ISABEL, CEM, INCL and LAQGSM. For 
all of the INC models, the following interaction stages are invoked: nuclear cascade, pre- 
equilibrium, evaporation or fission, and de-excitation. The user can choose to deactivate 
specific interactions. For results presented in this paper, the default option was used 
which invokes all interactions except for the pion off cases. For projectile energies E < 
3.5 GeV, the hadron physics modes Bertini, FLUKA, CEM03, INCL4 and LAQGSM are 
used. For heavy ion transport, LAQGSM is used exclusively. 

The MCNPX 2.6F physics package includes two significant changes. The recently 
improved version of the light ion event generator code CEM03, and the full implemen- 
tation of the LANL version of the LAQGSM, version 3.01 [26, 27, 28]. The CEM03 
model describes the nucleon induced reactions at incident energies up to 3.5 GeV [24, 29]. 
LAQGSM describes both the particle and nucleus induced reactions at energies above 5 
MeV and up to 1 TeV/nucleon for all isotopes. These two new improvements are im- 
portant to this study. LAQGSM impacts the ability to transport heavy ions and the 
improvements in CEM03 impact the light ion and pion transport. Detailed description of 
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Table 5: Physics model energy limits for heavy ions, light ions and hadrons. 



Heavy 

Light 

Hadrons 


ions 

ions 


1 

LAQGSM 

ISABEL (E <1 GeV) 

Bertini (E <3.5 GeV) 



LAQGSM ( E >1 GeV) 

FLUKA (E >3.5 GeV) 

2 

LAQGSM 

ISABEL (E < 1 GeV) 

CEM03 (E <3.5 GeV) 



LAQGSM (E >1 GeV) 

FLUKA (E >3.5 GeV) 

3 

LAQGSM 

INCL4 (E < 1 GeV) 

INCL4(£ <3.5 GeV) 



LAQGSM (E >1 GeV) 

LAQGSM {E >3.5 GeV) 

4 

LAQGSM 

LAQGSM 

Bertini (E <3.5 GeV) 
LAQGSM (E >3.5 GeV) 

5 

LAQGSM 

LAQGSM 

CEM03 (E <3.5 GeV) 
LAQGSM {E >3.5 GeV) 


the improvements in the LAQGSM and CEM03 event generation codes used in MCNPX 
2.6F can be found in Mashnik et al. [27, 28]. For the results described in this paper, the 
CEM03 model was used exclusively for all the hadrons physics with E < 3.5 GeV, and 
mix-match was used for E > 3.5 GeV. The option to force LAQGSM was tried; however, 
the code showed instability due to errors in handling some specific ion interactions. These 
instabilities were avoided by using the mix-match option. As shown in Table 5, using the 
mix-match option, either LAQGSM or FLLTKA is used for hadron energies E > 3.5 GeV. 
The mix-match algorithm dictates that LAQGSM physics is used first, and if it fails then 
it is passed on to FLUKA. Work is already underway to evaluate the use of LAQGSM 
exclusively, once these errors in LAQGSM are fixed and available for beta users. However, 
the exclusive use of LAQGSM is not expected to change the most significant findings de- 
tailed in this paper, because the GCR spectrum peaks below the transition energy for the 
LAQGSM physics model ( E > 3.5 GeV). 

3.3 Geometry and source 

The main goal of this paper is to quantify the impact of neutral and charged pion 
production and transport contribution on the downstream dose in human tissue for the 
1977 solar minimum GCR spectrum. Results for two different shield thickness of Al (5 
and 20 g/cm 2 ), followed by 30 g/cm 2 of human tissue, are discussed in this paper. These 
two cases are chosen because these are two standard benchmarks routinely used for V&V 
activities. We present particle flux and dose results, and compare the differences in these 
values with and without pion physics at various depths of Al and human tissue. The soft 
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tissue elemental composition was used from International Commission on Radiation Units 
(ICRU) 33 report [32], This section describes the details of the MCNPX input. 

A simple slab using a combinatorial geometry (CG) description is created for the two 
benchmark cases of 5 and 20 g/cm 2 of A1 shield followed by 30 g/crn 2 of tissue. Both 
cases are modeled individually to avoid any particle backscatter contributions and also to 
verify the results separately. The GCR source is modeled as a planar source uniformly 
distributed on a surface. It is modeled as an isotropic source, directionally biased towards 
the shield. Figure 1 shows the simulated geometry with the A1 and tissue thickness, and 
a wire frame GCR source plane centered on the left hand side. The dose segments are 
shown as solid colors, and the A1 and tissue slab is shown as transparent solids. 

The source is modeled explicitly with 26 heavy ions, helium (alpha) and hydrogen 
(proton) ions with the maximum kinetic energy of 5 GeV per nucleon. The source sam- 
pling rate is determined by the relative abundance of each element in the 1977 solar 
minimum GCR event. The GCR source boundary condition is calculated using the 2004 
Badhwar-O’Niell GCR model referred to as the B-O’04 Model [33]. Figure 2 shows the 
ion abundance curve for the 1977 solar minimum GCR event as predicted by the B-O’04 
model. Figure 3 shows the differential flux of the nine ions for the 1977 solar minimum 
GCR spectrum. 

3.4 Tallies 

MCNPX offers seven standard current and flux tally options. These, in combination 
with tally modifiers and multipliers, can be used to create a large selection of tally options. 
For the purpose of this analysis, a combination of four different tallies were used: surface 
current tally, flux tally, energy deposited per track length tally and collision heating tally. 
The cosine and energy multiplier cards and the “res” card were used to generate the 
results from MCNPX. The energy deposited per unit weight values were tallied at various 
A1 and tissue depths. These values are related to absorbed dose in units of cGy. The 
following conversion factor is used: 1 MeV per gram is equal to 1.6 x 1CU 8 cGy. For 
convenience, the energy per gram will be referred to as dose. 

4 Results and discussion 

All calculations presented in this paper are performed with MCNPX version 2.6F. 
For hadron transport E < 3.5 GeV, the CEM03 model is used exclusively; for kinetic 
energies E > 3.5 GeV, the mix-match option is used. For light ions, the mix-match 
physics option is used. The heavy ion transport is performed with the LAQGSM nuclear 
cascade/evaporation model. The calculations are performed to meet a minimum of 10% 
or less statistical uncertainty in the evaluated data. For all results, the error bars are 
plotted on the curves. In most cases, they are not visible because the relative error is very 
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small. The calculations are performed on a 60-node (120 CPU) cluster, “Poe”, at NASA 
Langley Research Center. Poe is a LINLIX cluster, with an AMD Dual core 2.4 GHz 
Opteron chip set with 2 GB dedicated memory per core. The MCNPX code is built with 
Intel 10.1 FORTRAN and a C++ compiler on a Fedora Core 4 Linux operating system. 

We will use the following notation for discussing these results. Two broad categories 
of results are presented in this paper: the flux and the dose results. For flux results, in 
most cases, the results are per source particle, rather than absolute values. Light ions 
(A < 4) include deuterons, tritons, helions and alphas. Heavy ions (A > 4) include all 
possible ions (2205 isotopes) heavier than alphas. Comparisons are made between flux 
and dose results for the pion channel turned on and off. The pion transport includes the 
charged and neutral pion channels. For all calculations, the neutron, photon, electron, 
light ion and heavy ion channels are included. The only difference between the pion on 
and pion off results is the presence or absence of the pion channel. 

First, we present results which show that the source and geometry were modeled 
correctly, and that the code behaved in the expected fashion with regard to the production 
of secondary particles and the attenuation of primaries. These are important checks to 
assure that differences observed in the final results are truly due to pion physics. We will 
then discuss the total particle flux as a function of depth in materials (A1 and tissue) 
and the particle differential fluxes. Finally, dose results are presented. These include 
dose-depth curves and relative dose contributions from various secondary products to the 
total dose. 

4.1 Source Verification 

Figures 4-6 show the source particle differential flux for proton, alpha and 16 0 at 0 
(boundary condition), 5, and 20 g/crn 2 of A1 depths. In the figures, the word “Source” 
refers to the boundary condition at 0 g/cm 2 . The zero depth flux curves show that the 
boundary condition is modeled correctly. These curves match the flux curves shown in 
Figure 3. All the source particles (28 in total) show similar agreement with the modeled 
boundary conditions. The flux curves for the 5 and 20 g/cm 2 of A1 shield show that the 
lower energy source particles attenuate relatively quickly as they pass through Ah The 
linear attenuation of charged particles show strong charge (Z) dependence. Hence, the 
attenuation is greater for 16 0 when compared to alphas and protons. Additionally, the 
lighter fragments (protons and alphas) are produced more frequently than heavier ions. 
Charged particle radiation of a specific energy has a maximum range in a given medium. 
At this range, all the particles are stopped regardless of source strength. At lower energies, 
we expect all the primary particles to be attenuated based on particle range. For example, 
the range of a 100 MeV proton is about 10 g/cm 2 in Al, whereas a 1000 MeV proton has 
a range of about 412 g/cm 2 in Al [34], The ranges for 100 and 1000 MeV alphas in Al 
are 0.85 g/cm 2 and 49 g/cm 2 , respectively. We see that proton flux above 1000 MeV is 
relatively unaltered when it passes through 5 and 20 g/cm 2 Al. The lower energy flux 
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is reduced considerably. Particles observed at lower energies at depths greater then the 
range of this energy particle are secondary particles. That is, the lower energy spectrum 
is a function of the production and loss of lower energy particles in the medium. Notice 
also that the error bars for the lower energy particles for alpha are large. The error bars 
for lower energy flux are larger for 16 0. This is because the frequency of production of 
these particles is lower which, in turn, impacts the convergence on these particles. 

4.2 Flux Results 

Figures 7-10 show the total particle flux versus depth curves for source particles 
(protons, alphas, heavy ions) and secondary light ions, neutrons, photons, and electrons. 
Figures 7 and 8 are flux versus depth curves in A1 and tissue for the case with pion 
channel on. Figures 9 and 10 are similar curves for the pion channel turned off. Source 
particle fluxes decrease as they attenuate through A1 and tissue. The flux attenuation 
rate is greater for heavy ions, as compared to alphas and protons. This is due to two 
reasons. First, the heavier particles attenuate more quickly due to Coulomb interactions. 
Second, a greater number of secondary protons and alphas are produced through the 
fragmentation process. The secondary particle build up is evident by observing flux 
versus depth curves for other light ions (deuterons, tritons and helions) which are not 
present in the primary source. We also note the build up of neutrons, photons, electrons 
and pions as a function of depth. We observe that the total flux for all the secondary 
products build up as the primary source particles travel through material. Note that 
these are not present in the source term, and hence are created entirely from the nuclear 
interactions and fragmentation of the primary particles. 

When comparing flux between A1 and tissues, we notice that the neutron flux builds up 
in Al, and it begins to reduce in tissue. This is expected since neutrons start to attenuate 
in tissue, due to the higher concentration of hydrogen in tissue composition. Note that 
as the neutrons attenuate, they produce light ion recoils. Notice the increase in deuteron 
and triton flux in tissue. 

The most interesting result from these plots is the difference in rate of build up of 
photon and electron fluxes in Al when compared between the pion on and pion off case. 
See Figures 7 and 9. Notice that for the pion off case, the photon curve curve lies on 
top of the neutron curve because most photons are produced from nuclear interactions. 
However, in the pion on case, the photon curve is much higher than the neutron flux. The 
increase in photon flux is due to the contributions from neutral pion decay, EM cascades 
and the charged pion interactions. These results relate to the theoretical discussion on 
mean lifetime and range of pions, as seen in Table 3. Notice that neutral pions have very 
short range and immediately decay into photons. Once produced, the charged pions travel 
deeper and undergo nuclear interactions. We notice in Figure 7 and 8 that photon and 
charged pion fluxes increase as a function of depth. The photon production from different 
interactions is clear when we examine the differential flux for neutrons and photons. 
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Figures 11 and 12 show the neutron and photon differential spectra. The neutron 
and photon plots are shown at the 5 and 20 g/cm 2 A1 depths. Additionally, neutron and 
photon spectra at 30 g/cm 2 tissue depth behind 20 g/cm 2 A1 are shown in the same plots. 
Notice that the neutron and photon spectral intensity is greater as the primary source 
particles pass through Al. The neutron and photon intensity for the 20 g/cm 2 case is 
greater than for the 5 g/cm 2 case. 

A close examination of Figure 12 reveals that there are three characteristic produc- 
tion regions in the photon spectra. The higher energy region ( E > 50 MeV), with a 
peak at about 100 MeV, is exclusively from neutral pion decay. The intermediate region 
(1 MeV - 50 MeV) is from nuclear interactions including neutron scattering reactions. 
Lastly, the low energy region (E < 1 MeV) is mainly from decay gammas and low energy 
Bremsstrahlung. Clearly, the pion produced photons, being energetic, would travel to 
significant depth. Consequently, these photons deliver radiation dose at greater depths. 
This is evident when we compare the photon spectra at 5 and 20 g/cm 2 of Al depths. 
The photon spectra remain practically unaltered at higher energies ( E > 100 MeV), and 
the photon flux increases in the intermediate energy range. 

Finally, let us examine the relative attenuation of neutrons and photons as they travel 
through 30 g/cm 2 tissue. Once in tissue, note that the lower energy neutrons (E < 
50 MeV) attenuate, but the high energy neutrons remain unaltered. Also, note that 
the 1 MeV neutron flux is reduced by an order of magnitude. However, in the case of 
photons, we notice that there is in fact a slight increase in high energy photons and only 
a small reduction in the lower energy photons. This significantly changes the make up of 
the secondary radiation products and their delivered dose as these photons travel through 
spacecraft and human tissue. Further characterization of the photon spectrum to evaluate 
its effect on dose equivalent is already underway. 

Tables 6 and 7 show the percent change in total particle flux due to the presence of 
pions in Al and tissue, respectively. It is clear from these tables that by including the pion 
transport, the particle flux for electrons and photons increase considerably, both in Al and 
in tissue. More importantly, as seen in Table 6, the total flux for neutrons changes steeply, 
(~ 7%) at 20 g/cm 2 in Al. The difference remains relatively flat (~ 6.5%) in tissue, as 
seen in Table 7. Similarly, the light ion flux increases sharply, (~ 3%) in Al and then 
remains flat throughout the depth of tissue. The proton flux increases slightly in both 
Al and in tissue. The heavy ion flux remain largely unaltered in Al and in tissue, when 
compared for the pion on and off cases. The change in total flux results are consistent 
with the theoretical discussion of the previous section. Due to pion interactions, we expect 
and notice an increase in particle fluxes for neutron, light ion and the electromagnetic 
particles (photons and electrons). Further study on kaons and muons is underway. 
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Table 6: Percent difference (between pion on and pion off) in total flux at various A1 
depths. 


Depth in 
Aluminum 
[cm] 

Electrons 

Photons 

Neutrons 

Ions (A < 4) 

Protons 

Ions (A > 4) 

1.0 

69.2% 

38.6% 

3.9% 

1.9% 

0.0% 

0.0% 

3.0 

74.5% 

41.4% 

4.7% 

2.3% 

0.1% 

0.0% 

5.0 

76.6% 

43.0% 

5.2% 

2.4% 

0.1% 

0.0% 

10.0 

80.9% 

46.1% 

6.1% 

2.3% 

0.2% 

0.0% 

15.0 

83.6% 

48.6% 

6.6% 

2.8% 

0.3% 

0.0% 

20.0 

86.0% 

51.0% 

7.0% 

2.7% 

0.4% 

-0.1% 


Table 7: Percent difference (between pion on and pion off) in total flux at various tissue 
depths. 


Depth in 
Tissue 
[cm] 

Electrons 

Photons 

Neutrons 

Ions (A < 4) 

Protons 

Ions (A > 4) 

1.0 

86.7% 

52.0% 

6.9% 

2.4% 

0.4% 

0.1% 

5.0 

87.2% 

54.7% 

6.6% 

2.2% 

0.5% 

0.1% 

10.0 

88.5% 

57.4% 

6.3% 

2.6% 

0.6% 

0.1% 

15.0 

89.7% 

59.8% 

6.3% 

2.3% 

0.7% 

0.3% 

20.0 

90.8% 

62.1% 

6.3% 

3.1% 

0.8% 

0.4% 

30.0 

93.1% 

66.5% 

6.2% 

2.7% 

1 .0% 

1.1% 


4.3 Dose Results 

Figures 13 and 14 show the relative dose contribution of various radiation constituents 
to the total dose at 30 g/cm 2 of tissue, after the primary radiation has passed through 
5 g/cm 2 A1 shielding for pions off and pions on, respectively. In the previous section, we 
noticed a change in flux due to the presence of pion physics. Let us now examine the 
effect on dose. Comparing Figures 13 and 14, we notice a significant difference in photon 
contribution to dose. It changes from about 1% for the pion off case, to about 9% for the 
pion on case, which is an increase by a factor of 9 (89%). The increase in dose is related 
to the high energy photons (E > 100 MeV) that are produced from the pion channel. 
The direct pion contribution (through charged pion interactions) is ~ 3%. The total dose 
contribution of pions and photons combined is ~ 13.5%. We do not see much change in 
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neutron and heavy ion contributions to the total dose. The dose contribution from light 
ion decreases by ~ 20%. For protons, the dose contribution decreases by ~ 10%. We 
note a change in makeup of the relative dose contributions to total dose from individual 
constituents. This is a significant and important finding. The implication of this change 
to dose equivalent is being investigated. 

Figures 15 and 16 show the results for the 20 g/cm 2 of Al. The change in dose trends 
for the 20 g/cm 2 of Al case are similar to the 5 g/cm 2 case. The changes are even greater 
in the 20 g/cm 2 of Al case. The photon dose contribution changes to 13% (pion on) 
from 1% (pion off). Hence, the photon dose increases by ~ 92%. The direct pion dose 
contribution is ~ 5%. The total combined contribution for photon and pion of ~ 18%, is 
significant. The heavy ion dose contribution remains the same between the pion on and off 
cases. The overall heavy ion dose contribution decreases from 3% to 2% when compared 
to the 5 g/cm 2 of Al case. This is expected because additional aluminium shielding breaks 
up and attenuates the primary heavy ions, resulting in a lower direct dose contribution 
from heavy ions at deeper depths. The relative dose contribution of neutrons to total 
dose changes from 5% to 4% (20% decrease). Similarly, the light ion dose contribution 
changes from 13% to 11% (15% decrease); whereas, the proton dose contribution changes 
from 79% to 65% (18% decrease). 

It is important to point out that the changes noted in the dose contributions are in the 
relative contribution of each constituent to the total dose. As we will see later, the total 
dose increases due to the pion channel, additionally, the dose from each of the constituents 
increases deeper into the tissue. Figures 13 - 16 showed us how the relative contributions 
change as the pion physics and transport are included. 

Figures 17 - 20 show the relative dose contribution to total dose results for the 5 g/cm 2 
and 20 g/cm 2 of Al at various tissue depths. These results show that there is clearly 
a significant change in the makeup of dose contributions to total dose from individual 
constituents, due to the presence of the pion channels. 

Figure 21 shows the percent change in dose, at various tissue depths, for each individual 
constituent due to the presence of the pion channels. The results shown are for 20 g/cm 2 
of Al. The trends for the 5 g/cm 2 of Al case are very similar and hence not shown here. 
Notice that the dose from light and heavy ions increases by about 5% at 30 g/cm 2 of tissue 
depth. The light ion dose increases by about 2% at the surface of the tissue (0 g/cm 2 ) 
with a gradual increase as a function of tissue depth. For heavy ions, the slope is bigger 
than for light ions from 0 g/cm 2 to 30 g/cm 2 tissue depth. The neutron dose changes by 
nearly ~ 9%, which is a noticeable increase. The slope for the neutron dose remains flat 
through the tissue depth, consistent with the neutron flux results. The contribution from 
protons changes by a smaller amount, ~ 2.5%. This relatively smaller change in proton 
dose can be explained. We note that the dose from protons is dominated by the primary 
source particles. Hence, the relative increase in proton dose due to the pion channel does 
not show a significant change, as it does for the neutrons. An important finding from these 
results is the overall increase in dose from hadrons and ions deeper into the tissue depth. 
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This is significant because by ignoring pions and EM particles, the dose contributions to 
critical organs could be under predicted. This could lead to misinterpretation of shielding 
effectiveness. Further investigation using linear energy transfer (LET) and dose equivalent 
is needed to fully evaluate the broader impact of these particles. 

Figure 22 shows the increase in total dose as a function of tissue depth behind 5 g/cm 2 
and 20 g/cm 2 Al. We observe that the increase in total dose is nearly 16% due to the 
presence of the pion, when the primary GCR spectrum travels through 20 g/cm 2 of Al and 
30 g/cm 2 of tissue. For the 5 g/cm 2 case, the total dose increases by ~ 9%. In summary, 
we notice a net increase in total dose when the pion channel is turned on. We observe 
that this increase is attributed to both increased nuclear interactions and EM cascades. 
An increase in dose as a function of tissue depth from all constituents is also noted. 


5 Conclusions 

Results of this study suggest that inclusion of meson, lepton and photon physics for 
space radiation shielding calculations may be necessary. Further investigation of this 
effect on dose equivalent and linear energy transfer (LET) is being pursued to quantify 
the effects of mesons on integral quantities. Both the flux and dose results suggest that 
inclusion of these particles will affect the overall dose prediction and the flux distribution 
of hadrons. This study shows a statistically significant difference in the contribution 
of absorbed dose, from the various secondary constituents, to the total absorbed dose 
when pion transport is included, as a substantial increase in photon flux and dose is 
observed. The photon spectrum shows three clear regions of particle production. A 
considerable portion of the spectrum is at energies greater than 100 MeV, which is a 
potential concern for shielding consideration. HZETRN has proved to be an efficient and 
a reliable tool for space radiation shielding evaluation for NASA. Since the underlying 
physics and transport are handled differently in the HZETRN and MC codes, specific code 
comparison benchmarks must be developed to investigate the impact of these findings. 
Currently, V&V efforts are underway to compare these results with other MC codes. 
These results suggest that the V&V efforts of the transport codes for space radiation 
must further evaluate the meson contribution to human and electronics exposure inside 
spacecraft. 
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Figure 1: Plot of MCNPX geometry. 



Figure 2: Galactic cosmic ray (GCR) ion abundance. 
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Figure 3: Differential flux of 1977 GCR ions. 
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Figure 4: Proton flux at boundary, 5 g/cm 2 and 20 g/cnr Al. 
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Figure 5: Alpha flux at boundary, 5 g/cm 2 and 20 g/crn 2 Al. 
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Figure 6: 16 0 flux at boundary, 5 g/crn 2 and 20 g/cm 2 Al. 
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Figure 7: Particle flux at various depths of A1 with pion channel on. 
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Figure 8: Particle flux at various tissue depths behind 20 g/cm 2 A1 with pion channel on. 
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Figure 9: Particle flux at various depths of A1 with pion channel off. 


-♦-neutrons -"-photons protons -e-deutrons hk- tritons 
-i-helions alphas -*-(A>4) electrons -&- pions 

1E+01 i 
A 

03 


7 1E+00 - 



IE-04 


0 5 10 15 20 

Depth in Tissue (g/cm 2 ) 

Figure 10: Particle flux at various tissue depths behind 20 g/cm 2 A1 with pion channel 
off. 
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Figure 11: Neutron flux behind 5 g/cm 2 , 20 g/cm 2 Al and at 30 g/cm 2 tissue. 
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Figure 12: Photon flux behind 5 g/cm 2 , 20 g/cm 2 Al and at 30 g/cm 2 tissue. 
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Figure 13: Relative contribution of various constituents to total dose behind 5 g/cm 
A1 and 30 g/cm 2 of tissue, with the pion channel off. 



Figure 14: Relative contribution of various constituents to total dose behind 5 g/cm 
A1 and 30 g/cm 2 of tissue, with the pion channel on. 
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Figure 15: Relative contribution of various constituents to total dose behind 20 g/cm 2 of 
A1 and 30 g/cm 2 of tissue, with the pion channel off. 
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Figure 16: Relative contribution of various constituents to total dose behind 20 g/cm 2 of 
A1 and 30 g/cm 2 of tissue, with the pion channel on. 
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Figure 17: Relative contribution of various constituents to total dose behind 5 g/cm 2 of 
Al, with the pion channel off. 
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Figure 18: Relative contribution of various constituents to total dose behind 5 g/cm 2 of 
Al, with the pion channel on. 
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Figure 19: Relative contribution of various constituents to total dose behind 20 g/cm 2 of 
Al, with the pion channel off. 
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Figure 20: Relative contribution of various constituents to total dose behind 20 g/cm 2 of 
Al, with the pion channel on. 
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Figure 21: Percent increase in dose deposited by radiation constituents at various depths 
of tissue behind 20 g/cm 2 of Al. 



Figure 22: Percent increase in dose deposited at various depths of tissue behind 5 g/cm 2 
and 20 g/cm 2 of Al. 
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